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Abstract 

In this work, the quantization of the most general Bianchi Type I geometry, 
with and without a cosmological constant, is considered. In the spirit of identifying 
and subsequently removing as many gauge degrees of freedom as possible, a reduc- 
tion of the initial 6-dimensional configuration space is presented. This reduction 
is achieved by imposing as additional conditions on the wave function, the quan- 
tum version of the -linear in momenta- classical integrals of motion (conditional 
symmetries). The vector fields inferred from these integrals induce, through their 
integral curves, motions in the configuration space which can be identified to the 
action of the automorphism group of Type I, i.e. GL(3, 5ft). Thus, a wave function 
depending on one degree of freedom, namely the determinant of the scale factor 
matrix, is found. 

A measure for constructing the Hilbert space is proposed. This measure respects 
the above mentioned symmetries, and is also invariant under the classical property 
of covariance under arbitrary scalings of the Hamiltonian (quadratic constraint). 
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1 Introduction 



Since the conception by Einstein of General Relativity Theory, a great many efforts have 
been devoted by many scientists to the construction of a consistent quantum theory of 
gravity. These efforts can de divided into two main approaches: 

(a) perturbative, in which one splits the metric into a background (kinematical) part 
and a dynamical one: = r/^ + h^ and tries to quantize h^. The only conclusive 
results existing, are that the theory thus obtained is highly nonrenormalizable [I]. 

(b) non perturbative, in which one tries to keep the twofold role of the metric (kinemat- 
ical and dynamical) intact. A hallmark in this direction is canonical quantization. 

In trying to implement this scheme for gravity, one faces the problem of quantizing a 
constrained system. The main steps one has to follow are: 

(i) define the basic operators "g^ and 7r M1/ and the canonical commutation relation 
they satisfy. 

(ii) define quantum operators whose classical counterparts are the constraint func- 
tions H^. 

(iii) define the quantum states ^[g] as the common null eigenvector of H^, i.e. those 
satisfying H^[g] = 0. (As a consequence, one has to check that H^, form a closed 
algebra under the basic Canonical Commutation Relations (CCR).) 

(iv) find the states and define the inner product in the space of these states. 

It is fair to say that the full program has not yet been carried out, although partial 
steps have been made 0. Concerning point (iii) we deem it pertinent to clarify the 
meaning of the imposition of the quantum constraints upon *S?[g]. A straightforward 
(modulo regularization prescriptions) but tedious calculation shows that any functional 
which is not a scalar functional of the curvature invariants (see ||) does not solve the 
linear constraints. Therefore, the imposition of the linear constraints, ensures that the 
wave functional will be a (scalar) functional of the 3-geometry and not of the coordinate 
system. Then, the dynamical evolution is provided by the quadratic constraint; the 
consistency of the quantum algebra, guarantees that the final wave functional, will be 
independent of the 4 dimensional coordinate system. 

In the absence of a full solution to the problem, people have turned to what is 
generally known as quantum cosmology. This is an approximation to quantum gravity 
in which one freezes out all but a finite number of degrees of freedom, and quantizes the 
rest. In this way one is left with a much more manageable problem that is essentially 
quantum mechanics with constraints. Over the years, many models have appeared in the 
literature In most of them, the minisuperspace is flat and the gravitational field is 
represented by no more than three degrees of freedom (generically the three scale factors 
of some anisotropic Bianchi Type model ||). 
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In a series of earlier publications ||, we have considered the quantization of the 
most general spatially homogeneous spacetime for various Bianchi Types. Thus, our -in 
principle- dynamical variables were the scale factors 7 a /j(i), the lapse function N(t) and 
the shift vector N a (t). 

The presence of the linear constraints -along with the conditional symmetries of the cor- 
responding Hamiltonian- enabled a reduction of the initial configuration space to a lower 
dimensional one, spanned by the curvature invariants characterizing the 3-geometry. The 
ultimate justification of this reduction is the fact that -from the point of view of the 
3-geometry- the omitted degrees of freedom, are not physical but gauge J7|. 

The case of Bianchi Type I geometries, has been repeatedly treated in the literature - 
both at the classical level || and the quantum level |§. In all these works, the scale factor 
matrix is taken to be diagonal. The role of the topology of the 3-slices, is emphasized in 
the work of Ashtekar and particularly in that of Hervik where a non trivial T 3 topology 
leads to a nine dimensional moduli space. The main reason for this plethora of treatments 
of Type I, is the simplicity brought by the vanishing structure constants, i.e. the high 
spatial symmetry of the model. It is true that at the classical level, the scale factor 
matrix, can be diagonalized on mass-shell -through a constant matrix [|Kj- while the 
shift can be set equal to zero. However, if one intends to give weight to all states, 
one has to start with the most general form which is described by the 6 scale factors 
lapif) an d the lapse function N{t). The absence of H a s due to the vanishing of the 
(7g 7 s, implies that in principle all 7 aj «s are candidates as arguments for the wave function 
which solves the quadratic constraint (Wheeler-DeWitt equation). This is in contrast to 
what happens in other Bianchi Types f§ where, less or equal to 3, combinations of 7q,/3S 
and C^ 7 s, parameterize the reduced configuration space. 

In this short communication, we present a complete reduction of the initial con- 
figuration space for Bianchi Type I geometry -by extracting as many gauge degrees of 
freedom, as possible. Two separate cases are considered; when the cosmological constant 
is present and when is not. In either case, a wave function which depends on one degree 
of freedom (namely the determinant of the scale factor matrix) is found, by imposing on 
it, the quantum versions of all classical integrals of motion as additional conditions. 

2 Bianchi Type I, the A ^ case 

In this work, we will quantize the known action corresponding to the most general Bianchi 
Type I cosmologies, i.e. the action giving Einstein's Field Equations derived from the 
line element: 

ds 2 = (-N 2 (t) + N a (t)N a (t))dt 2 + 2N a (t)a?(x)dx i dt + -f al3 (t)a^(x)a^(x)dx i dx j (2.1) 

where a" 1 are the invariant basis one-forms of the homogeneous surfaces of simultaneity 
S 4 , satisfying: 

do" = C a M a?Ao^ a° j -o? ti = 1C% a] of (2.2) 
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with C^ 7 being the structure constants of the corresponding isometry group. In 3 di- 
mensions, the tensor admits a unique decomposition in terms of a contravariant 
symmetric tensor density of weight —1, m al3 , and a covariant vector u a = as fol- 

lows: 

C% = m aS es Pl + Vf,6° - u 7 5%. (2.3) 
For the Class A {y a = 0) Bianchi Type I, this matrix is [II : 




m af3 = I | (2.4) 

resulting in vanishing structure constants. 

As is well known [12], the Hamiltonian of the above system is H = N(t)H +N a (t)H a , 
where: 

H = \L a(3 ^^ + 1 R + 1 K (2.5) 



is the quadratic constraint, with: 



Laf3fj,u la/ilfti/ lavlflfi lafllfiv 



(2.6) 



7 being the determinant of i a/ 3, and: 

Ha = C^r P y (2.7) 

are the linear constraints. Note that N appearing in the Hamiltonian, is to be identified 
with N/ y^y. For all Class A Types, the canonical equations of motion, following from 



75]), are equivalent to Einstein's equations derived from line element (|2.1| ) -see [|l2jj . 
The quantities H , H a are weakly vanishing [HJ, i.e. H « 0, H a « 0. For all Class 
A Bianchi Types (C^g = 0), it can be seen -using the basic PBR {i a p, ir^} = 5^- that 
these constraints are first class, obeying the following algebra 

{Ho, H a } = ^ ^ ^ 

{H a ,Hp} = —^CZpHy, 

which ensures their preservation in time, i.e. Hq ~ 0, H a « 0, and establishes the 
consistency of the action. 

If we follow Dirac's general proposal |13[ for quantizing this action, we have to turn 
H , H a , into operators annihilating the wave function 

In the Schrodinger representation: 

la/3 -> lafi = la/3 , , 

n a0 ^ -a/3 = _ % 8 (2.9) 
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with the relevant operators satisfying the basic Canonical Commutation Relations (CCR) 
-corresponding to the classical ones: 

i%^n=<e=~m+s^a)- (2.10) 

In Bianchi Type I (Cp = 0), the second of Q2.6|) gives R = 0, while relations (|2] 
vanish identically, and the algebra (|2.8| ) is trivially satisfied. Thus, ( |2.5|) , reads: 



ffo-V^ + TA (2.11) 

Thus, the only operator which must annihilate the wave function, is Hq; and the 
Wheeler-DeWitt equation = 0, will produce a wave function, initially residing on a 
6-dimensional configuration space -spanned by s. The discussion however, does not 
end here. If the linear constraints existed, a first reduction of the initial configuration 



space, would take place [14|. New variables, instead of the 6 scale factors, would emerge 



-say q l } with % < 6. Then a new "physical" metric would be induced: 

g ij = L a ^ v dql dq3 (2.12) 



According to Kuchaf's and Hajicek's [14[] prescription, the "kinetic" part of Hq would 



have to be realized as the conformal Laplacian (in order for the equation to respect the 
conformal covariance of the classical action), based on the physical metric ( |2.12| ). In 
the presence of conditional symmetries, further reduction can take place, a new physical 
metric would then be defined similarly, and the above mentioned prescription, would 
have to be used after the final reduction || |15| . 

The case of Bianchi Type I, is an extreme example in which all the linear constraints, 
vanish identically; thus no initial physical metric, exists -another peculiarity reflecting 
the high spatial symmetry of the model under consideration. In compensation, a lot 
of integrals of motion exist ant the problem of reduction, finds its solution through the 
notion of "Conditional Symmetries". These linear in momenta integrals of motion, if 
seen as vector field on the configuration space spanned by 7 Q /3S, induce -through their 
integral curves- motions of the form 7 Q/3 = A^A^^, A e GL(3, 3ft) (see section 2 of 0) 
which not only are identical to the action of spatial diffeomorphisms, but also describe 
the action of the automorphism group -since GL(3, 3ft) is the Aut(G) for Type I |16| . 

The generators of this automorphism group, are (in a collective form and matrix 
notation) the following 9 -one for each parameter: 




\fo= | e C V | , IE [1,... ,9] (2.13) 
with the defining property: 

C^K = C^K + C^. (2.14) 
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Exponentiating all these matrices, one obtains the outer automorphism group of Type 
I, since there is not Inner Automorphism subgroup (all structure constants vanish). 

For full pure gravity, Kuchaf |L5] has shown that there are no other first-class func- 
tions, homogeneous and linear in the momenta, except the linear constraints. If however, 
we impose extra symmetries (e.g. the Bianchi Type I -here considered), such quantities 
may emerge -as it will be shown. We are therefore -according to Dirac [0- justified 



to seek the generators of these extra symmetries; their quantum-operator analogues will 
be imposed as additional conditions on the wave function. The justification for such an 
action, is obvious since these generators correspond to spatial diffeomorphisms -which 
are the covariance of the theory. Thus, these additional conditions are expected to lead 
us to the final reduction, by revealing the true degrees of freedom. Such quantities are, 
generally, called in the literature "Conditional Symmetries" [15]. 






From matrices ( |2.13|) , we can construct the linear -in momenta- quantities: 

Em = Xfofffap** (2.15) 
In order to write analytically these quantities, the following base is chosen: 

Ai = I I , A 2 = 

A 5 = I I , A 6 = ( 1 ] (2.16) 
\l o o) 
1 \ 
A 7 = I -1 I , A 8 = -1 , A 
\ / \ 1 / 

It is straightforward to calculate the Poisson Brackets between Em and H : 

{E (Ih H } = - 1 AX a a (2.17) 

But, it holds that: 

E m = {E (I) ,H } = - 1 A\ a a (2.18) 

-the last equality emerging by virtue of ( |2.17| ). Thus: 

Em = {Em, H } = => Em = Km = constants, I e [1, ... ,8] (2.19) 

We therefore conclude that, the first eight quantities Em, are first-class, and thus inte- 
grals of motion. Out of the eight quantities Em, only five are functionally independent 
(i.e. linearly independent, if we allow for the coefficients of the linear combination, to 
be functions of the j a p s); numerically, they are all independent. 
The algebra of Em can be easily seen to be: 

{Em,E {J) } = ~C?jE {M) , I, J, Me [1,... ,9] (2.20) 
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where: 



(J)\ 



Wj A (M), 



y,Me[i,...,9] 



(2.21) 



the square brackets denoting matrix commutation. 

The non vanishing structure constants of the algebra ( |2.21| ), are found to be: 



°13 


= 1 


°15 — 


-1 


°16 


= 1 


r l - 

°18 - 


-1 


r 1 

°17 


= -2 


°24 


= 1 


°25 — 


1 


°25 


= -1 


^26 = 


-1 


^2 
°27 


= -1 




= 1 


°34 — 


-1 


C35 


= 1 


C37 = 


1 


°38 


= 2 




= 1 


°47 — 


-1 


^4 
°48 


= -2 


^57 = 


1 


°58 


= -1 


^67 


= 2 


°68 — 


1 















(2.22) 



At this point, in order to achieve the desired reduction, we propose that the quantities 
Em -with / G [1, . . . , 8]- must be promoted to operational conditions acting on the 
requested wave function \l/ -since they are first class quantities and thus integrals of 
motion (see ( |2.19| )). In the Schrodinger representation: 



a/3 



IE [1,... 



(2.23) 



In general, systems of equations of this type, must satisfy consistency conditions decreed 
by the Frobenius Theorem: 



K {J) K {I) ^ 



Subtraction of these two and usage of Q2.20| ), results in: 

=> C¥ T K(M\ = 



(2.24) 



(2.25) 



This relation constitutes a selection rule for the numerical values of the integrals of 



motion. In view of the Lie Algebra ( |2.22| ), selection rule (|2.25|) imposes that all Ks 
vanish, i.e. K\ — . . . — Kg = 0. This fact restores the action of the diffeomorphisms as 
covariances of the quantum theory, in the sense that now, we have conditions of the form 
E(i}ty = 0. Instead, if we also had E<§\ (as is the case A = 0) then K 9 would remain 



arbitrary. With this outcome, and using the method of characteristics flT7|, the system 



of the five functionally independent P.D.E. s ( |2.23|) , can be integrated. The result is: 

*(7) 



(2.26) 



i.e. an arbitrary (but well behaved) function of 7 -the determinant of the scale factor 
matrix. 

A note is pertinent here; from basic abstract algebra, is well known that the basis of 
a linear vector space, is unique -modulo linear mixtures. Thus, although the form of the 



system ( \2. 23\ ) is base dependent, its solution ( \2.26] ), is base independent 
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The next step, is to construct the Wheeler-DeWitt equation which is to be solved by 
the wave function ( |2.26| ). The degree of freedom, is 1; the q = 7. According to Kuchaf's 



proposal fl4| , upon quantization, the kinetic part of Hamiltonian (|2.11|) is to be realized 
as the conformal Beltrami operator - based on the induced physical metric -according 
to ( gT2[) , with q = 7: 



11 r d l d l r JlnfS^v first of ( 2 - 6 ) o ,2 /« 97 x 



In the Schrodinger representation: 



- (2.28) 



where: 



□ 2 = D 2 = — 9 7 {V^I S n <9 7 } (2.29) 

is the 1-dimensional Laplacian based on gu (g ll gn = 1)- Note that in 1-dimension the 
conformal group is totally contained in the G.C.T. group, in the sense that any conformal 
transformation of the metric can not produce any change in the -trivial- geometry and 
is thus reachable by some G.C.T. Therefore, no extra term in needed in ( |2.29| ), as it can 



also formally be seen by taking the limit d — 1, R = in the general definition: 



Thus: 



H ^H = -I(-3 7 2 |5- - 3 7 |-) + A 7 (2.30) 



So, the Wheeler-DeWitt equation -by virtue of ( |2.26|) -, reads: 



H q = 7 V + + - 7 A^ = (2.31) 
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The general solution to this equation, is: 



*( 7 ) = Cl J (2^^) + c 2 Fo(2^/^) (2.32) 

where J n and Y n , are the Bessel Functions of the first and second kind respectively -both 
with zero argument- and Ci,C2, arbitrary constants. 

Some comments on this wave function. Indeed, at first sight, the fact that \1/ depends 
only on one argument and particularly on 7, seems to point to some undesirable degen- 
eracy regarding anisotropy; classically 7 can be gauged to e* and thus it seems as though 
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the anisotropy parameter does not enter \l/ at all. If, however, we reflect thoroughly, we 
will realize that this objection rests strongly on a -not generally accepted- mingling of 
the classical notion of anisotropy and the interpretation of the wave function. Indeed if 
we adopt the interpretation that the wave function \l/ (along with a suitable measure), is 
to give weight to all configurations parameterized by 7^, then the anisotropic configura- 
tion will, in general, acquire different probabilities. The degeneracy occurs only between 
two different anisotropic configurations with the same determinant 7. In compensation 
the scheme proposed here, avoids the gauge degrees of freedom as much as possible. The 
final probabilistic interpretation must await the selection of a proper measure, a theme 
which will be discussed in section 4. 



3 Bianchi Type I, the A = case 

This section mimics the previous one, with minor changes. Indeed, if the cosmological 
constant A is zero, some changes will take place. 

The first concerns the obvious alteration to the Hamiltonian ( [2 . 1 1|) , which now reads: 



H = -L af3fiu 7i a ^^ (3.1) 



This consequently, causes an alteration to the Poisson Bracket ( [2.17| ), which takes the 
form: 

{E {I) ,H } = } /€[!,..., 9] (3.2) 



while, ( |2.20|) still holds. Thus in the present case, there are nine integrals of motion 



-instead of eight. Also, the P.D.E. system ( |2.23| ) consists of nine members (instead of 



eight), but now out of the nine quantities Ey), only six are functionally independent; the 
previous five, plus the £V 9 ). Again, using the method of characteristics |T7|, the system 



of the six functionally independent P.D.E. s fl2.23|) , can be integrated. The result is: 



q = ci7^ 9/3 (3.3) 

where 7 is the determinant of the scale factor and Kg, the remaining constant -according 
to selection rule ( |2.25| ). 



The fact that this wave function does not depend on any combination of 7 Q/ 3S in 
an arbitrary manner (i.e. \l/ is not an arbitrary function of 7 a/ gs), might be taken as 
an indication that no reduced Wheeler-DeWitt equation can be written. On the other 
hand, this wave function does contain an arbitrary (essential) constant, which ought 
to be fixed by the dynamics. The puzzle can be solved by the following compromise; 
the initial configuration space, should be the mini-superspace i.e. we should write the 
Wheeler-DeWitt equation, based on the supermetric L a/3fJ-l/ . 

In the Schrodinger representation: 

l -L a ^^ -> ~nl (3.4) 



9 



Thus using ( |A.9| ), ( |A.11| ), ( |A.12| ) for d = 3 and V — 6, one may find respectively -see 
appendix: 



R 



15 



(3.5) 



L al3llu T^ u = -37 kA (3.6) 

and: 

9 , (9 (9 3 . i 

□c = L <*p^ + 37 KA ^ + 3.7 

Then Kuchaf's proposal for (|3.1|) reads: 



- 1 / 9 <9 3" 
H ^ H = -- [ L a/3fM/ - h 37 kA - h - ) (3.1 

2 V 0-f a fS-f^ #7 k a 2/ 



Substitution of the wave function Q3.3|) in the Wheeler-DeWitt equation Hq^/ = 0, 
with i^o given by the previous relation, determines the constant Kg. The outcome is: 

* = Cl7 V2/2 + (3. 9 ) 

The constants ci,C2, remain arbitrary and may be fixed after the selection of a proper 
measure via normalizability requirements. 



4 Measure and Probabilistic Interpretation 

In general, the issue of the selection of a measure, is an open question. For us, an 
important element for selecting it, is the conformal covariance; the supermetric L af3fiu is 
known only up to rescalings, because instead of N(t) one can take any N(t) = N(t)e~ 2uJ 

(with uj = Lo(^ a p)) and consequently L af3fiu (t) = L al3 ^ u (t)e 2ul . This property, is also 
inherited to the physical metric ( [2.12|) and is the reason for the Kuchaf's recipe, adopted 
in this work. 

It is therefore mandatory for the proposed measure /x, to be such that the probability 
density /i|\I/| 2 , be invariant under these scalings. Recalling that \I> = \I/e( 2_:D ) W//2 ( we 
conclude that fx must scale as ~p = fjLe^ v ~ 2 ^ u . 

It is not difficult to imagine such a quantity: any product of E(j) a p with E(j^ u (where 
E(i) a f3 = l/2(A^7 Kj a + (a <-> (3)) are the components of the vector field inferred from the 
integrals of motion £"(/)) has the desired property. Indeed the E^ s do not scale at all, 
while the supermetric scales as mentioned before. The group metric 0/j = Cf s Cj F can 
serve to close the group indices of Em a p. So, we arrive at the quantity: 

£ = L a ^e IJ E {I)a/3 E {J)lMU (4.1) 



10 



(where /J is the inverse of the group metric). The quantity £ transforms as: £ = £e 2 ^. 
Thus we only need to take jx = £( :D_2 )/ 2 Using the Lie algebra (|2.22| ), one obtains: 



©ij = Cf s Cj F 



/000006 
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6 
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12 
6 
0/ 



(4.2) 



Thus, in order to have an inverse, we restrict ourselves to the non trivial 8x8 subspace, 
and have: 
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) 



(4.3) 



That means that in the sum of the expression ( ^4.1| ), we exclude the vector Eg. After a 
straightforward calculation, one finds that: 



5 

12 



(4.4) 



The quantity £( I,_2 )/ 2 in each of the two cases (i.e. with and without the cosmological 
constant), defines the final expression for the measure fi. 

A drawback of the current measure, is the loss of the hermiticity of the Wheeler- 
DeWitt operator. This however is not such a blander because, the eigenvalues of this 
operator, are zero and thus loss of hermiticity does not result in having complex eigen- 
values. On the other hand the following welcomed features are obtained: 



In the case A ^ 0, the wave function which results in normalizability, is: 

= Ci J (2-y/^!p), where J is the Bessel function of the first kind with zero 
argument. Since this function has a branch cut discontinuity in the complex z 
plane running from — oo to 0, we are naturally led to consider A to be positive 
definite -since 7 is positive definite. Also, from the behaviour of the probability 
density, it is inferred that models with small A are much more probable. The limit 
A — > 0, gives a constant wave function and therefore constant probability. 



11 



• In the case A = 0, the wave function turns out to be non normalizable -if we use the 
measure above. This can be thought as an indication against the favour of such a 
model over the previous one. Of course, one could conceive of measures, e.g. gaus- 
sian ones, which would make the probability normalizable. Further investigation 
on this issue is irrelevant in view of the plethora of these possibilities. 

5 Discussion 

The present work is the last in a series of publications in which stress has been laid 
on the spirit and the usage of "Conditional Symmetries" within the programm of the 
Quantum Cosmology approximation towards the quantization of Gravity. It was shown 
how these can be used in order to achieve a reduction of the initial configuration space 
-a problem closely related to the symmetry of a model in general, and in particularly 
to the high spatial symmetry of Bianchi Type I model. Indeed, when these Conditional 
Symmetries are imposed on the requested wave function, i.e. the P.D.E. systems (|2.23| ), 
along with their consistency conditions, a significant reduction of the configuration space 
takes place, by extracting as much gauge degrees of freedom, as possible. This is a sine 
qua non procedure in view of the fact that these conditional symmetries in Quantum 
Cosmology, are nothing but a part of the G.C.T. group as is explicitly shown in 0. At 
the same time, most of the classical constants of motion, are set equal to zero -a fact 
not encountered before. 

The reduced space, is one dimensional and the degree of freedom left is still gauge (7 is a 
density) if seen from the 3-slice point of view. In both of the cases under consideration, 
a wave function is gained and a kind of "proposed role" for it, is exhibited. Only in 
the first case, the wave function is normalizable. Also, in this case, the results are in 
agreement with the latest estimations, concerning A, i.e. that it is positive and relatively 
small. Another peculiarity is that in A 7^ case, there is a non zero probability when 
7^0 but this it might be seen as a side-effect of the "preferred" coordinate system. 
Indeed, when 7 — > 0, the same limit is attained by the determinant of the supermetric, as 
well. If however, one forms any metric invariant in the minisuperspace, he will see that 
the initially detected singularity, is apparent and is due to the choice of the coordinate 
system of this space, only. 

At this point a remark concerning the implicit functional form of wave functions with 
respect to the 3-metric, is pertinent. In the literature it is common to assume a wave 
function dependant upon a diagonal metric. Such an assumption is certainly under- 
standable from the a classical point of view, but not justified quantum mechanically for 
two reasons: 

Ri The wave function is supposed to give weight to all states, not only to the classical 
ones. So, ab initio, it must depend on all 6 scale factors y a p. 

R2 As it has been mentioned in [|7|, plj, there is a particular class of G.C.T.s, which 
preserve manifest homogeneity of the line element of the generic Bianchi Homo- 
geneous 3-space, has a non trivial action on the configuration space spanned by 
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7a / s' s ) an d this action is that of the Automorphism group. Its differential descrip- 
tion, leads to the vector fields ( |2.15| ). Classically, in some Bianchi Types (e.g. Type 
VIII, IX) an off mass-shell diagonalization of the scale factor matrix is possible, 
and in some others (e.g. Type /) the diagonalization is possible only on mass-shell 
[[Kill . If one passes to the quantum level, i.e. if one uses the generators A^ of the 
corresponding simplifying automorphisms A", to construct vector fields of the type 
^plao-jpf- , he will encounter a very different situation; imposing these fields upon 
\I> depending on all 7 a /?s, he will not get a \& depending on the {711, 722, 733} only. 
For example the automorphism group for Bianchi Type IX is 5*0(3), so A" are 
the rotations, and their generators, are: 



Ai = -1 , A 










-n 
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1 
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-1 


/ 



The corresponding vector fields ( |2.15| ) upon ^ do not imply that this in prin- 
ciple, an arbitrary function of {711,722,733} only -while at the classical level, the 
arbitrary scale factor matrix is diagonalized through the usage of three independent 
rotations. 

On the issue of the measure, we have succeeded in constructing a measure which respects 
the scalings of H and gives a square integrable wave function (in the A 7^ case), at 
the cost of loosing hermiticity of the Wheeler-DeWitt operator. The last is not a serious 
drawback of the whole scheme, because not only the Wheeler-DeWitt equation, consti- 
tutes a zero eigenfunction problem, but also, the eigenfunctions (i.e. the permissible 
wave functions) are complex functions in general -the probability density in defined as 
the measure times the norm of the wave function. 

Lastly we would like to point out that the usage of the above techniques, is in contact 



with the work \TE\ where the strong gravity limit is treated. The connection to the strong 
gravity limit of all the other Bianchi Types, can be made much more explicit. The key 
observation is that the automorphism groups of all other Bianchi Types, are subsets of 
GL(3, 3ft) i.e. the automorphism group of Bianchi Type I. Consequently, if we adopt the 
point of view of imposing the appropriate subset of the linear integrals (|2.15|) on the wave 
function, corresponding to the automorphism group of the desired Bianchi Type, we will 
get a wave function which will satisfy the strong gravity limit of the Wheeler-DeWitt 
equation. For example, if we select the subset (in collective form) with: 

(k + /i x y \ 
k p ) (5.1) 
a n ) 

we will get a \l/ depending on 7 and q = 7^/27 and thus arrive at the strong gravity 

limit of equation (19) in the first of ||. Likewise for all Bianchi Types. 

This line of reasoning may also satisfy those which, driven by anisotropy "classical" 
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considerations, would like to have more than one arguments in the wave function. Of 
course, our point of view is that the less gauge freedom in the wave function the better 
and hence we adopt the procedure all the operator analogues of the integrals of motion 
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A Appendix 

In this appendix, we give some useful formulae, concerning the mini-superspace. 

Using the results of canonical analysis in a (d + l)-dimensional manifold, endowed 
with the line element ( |2.1|) , one arrives at the notion of mini-superspace spanned by 7 a/ g 
s (co-ordinates), and having as covariant metric the following: 

jj*Pv» = I( 7 ^ 7 ^ + 7 Q ^ 7 ^ _ 27 a/3 7^) (A.l) 
while the contravariant metric, is defined as: 

2 

L aPlw = {laplf3u + lavlflp ~ J—^l a plpv) (A. 2) 



in the sense that: 



L a ^ x L KXflu = V* = l - (Sffi + 5 a X) (A-3) 



The Christoffel symbols are defined as: 



T a f/ V = -L KXpa {L pa ^' al3 + L al3pa '^ u - L a ^' pa } (A.4) 



where: 

ajappv 

jjy.fi pv, pa _ 



Combined usage of (|A~l| ), (|A]|), ( [OP and ([Aj), gives: 



In the same spirit, the Riemann tensor is defined as follows: 

^af3pa-pu)^pu -pctf3 
u>£ kX w£ 



•QCtftpupv -pa{3pa,pv -pct0pi',pa , -paf3pa-puj^pu -paftpiz-puj^pcr / \ r>\ 

U k\ — 1 kA _ 1 kA + 1 w£ 1 kA ~ 1 w£ 1 kA 



where: 

-paf3pv,pa _ UL K \ 

kA " d lpa 
Contraction of (p, a) with (ft, A) results in the Ricci tensor: 

— L kX kA X K A J-kA l A -'J 
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A lengthy but straightforward calculation, gives: 

= hd^^ + dj au ^ - 2 7 Q V") (A.8) 
8 

With the help of ( |A.8|) and (|A.2| ) the Ricci scalar is found to be: 



R = L a ^ u R a/3fiV = ^(d 3 + d 2 - 2d) (A.9) 
Finally, the conformal Beltrami opetaror, is: 

where V is the dimension of the general metric space: V = , i.e. the number of the 
independent 7^. 
One can find that: 

W^:r~7a (All) 



thus ( |A.10| ), takes the form: 

2 T d 2 3-d 2 d V-2 „ 

frlafilp, rf -l #7«a 4(£>-l) 
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